Ab initio study of magnesium alanate, Mg(AlH 4 ) 2 
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Magnesium alanate Mg(AlH 4 ) 2 has recently raised interest as a potential material for hydro- 
gen storage. We apply ab initio calculations to characterize structural, electronic and energetic 
properties of Mg(AlH 4 ) 2 . Density functional theory calculations within the generalized gradient 
approximation (GGA) are used to optimize the geometry and obtain the electronic structure. The 
latter is also studied by quasi-particle calculations at the GW level. Mg(AlH 4 ) 2 is a large band 
gap insulator with a fundamental band gap of 6.5 eV. The hydrogen atoms are bonded in AIH4 
complexes, whose states dominate both the valence and the conduction bands. On the basis of total 
energies, the formation enthalpy of Mg(AlH 4 ) 2 with respect to bulk magnesium, bulk aluminum 
and hydrogen gas is 0.17 eV/Hb (at T = 0). Including corrections due to the zero point vibrations 
of the hydrogen atoms this number decreases to 0.10 eV/Hb. The enthalpy of the dehydrogenation 
reaction Mg(AlH 4 ) 2 — > MgH 2 + 2A1 + 3H2 (<?) is close to zero, which impairs the potential usefulness 
of magnesium alanate as a hydrogen storage material. 

PACS numbers: 61.50.Lt, 61.66.Fn, 71.20.Nr 
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The interest in alanate compounds as hydrogen stor- 
age materials was recently rekindled as the kinetics of 
hydrogen adsorption/desorption improved dramatically 
by the addition of transition metal catalystsi^ Alanates 
MAIH4, with M a lightweight alkaline metal, have a high 
gravimetric hydrogen density, which is essential for their 
application as storage materials. Most attention up to 
now has gone to sodium alanate, NaAlfi4, which has a 
hydrogen capacity of 7.5 wt. It releases hydrogen 

in a two stage process. The two stages involve reaction 
enthalpies that are sufficiently small to be of interest, 
namely 0.38 and 0.34 eV per H 2 molecule respectively. 
However, only three out of four hydrogen atoms are re- 
leased in this process, which lowers the effective hydrogen 
storage capacity. This has stimulated the search for other 
suitable alanates. 

Alanates M(A1H 4 ) 9 , with M a lightweight alkaline 
earth metal, have an even higher gravimetric hydrogen 
density. Recent interest turned to magnesium alanate, 
Mg(AlH 4 ) 2 , which has a hydrogen capacity of 9.3 wt. 
%i 5 i 6 i 7 i 8 Upon heating Mg(AlH 4 ) 2 releases hydrogen ac- 
cording to the reaction 



Mg(AlH 4 ) 2 -> MgH 2 + 2A1 + 3H 2 (<?). 



(1) 



Since decomposition of MgH 2 takes place at too high a 
temperature to be of practical usef2iS it is the amount of 
hydrogen released in QJ that determines the actual stor- 
age capacity of Mg(AlH 4 ) 2 . Still, this relatively large 
amount of 7.0 wt. % makes magnesium alanate a good 
candidate for hydrogen storage. Only little is known 
about the thermodynamics of this material, however^S 
Since up to now its synthesis has proceeded via an in- 
direct route, the first question is whether Mg(AlH 4 ) 2 is plane wave basis set. The reaction enthalpies are calcu- 



thermo dynamically stable with respect to decomposition 
into its elements. The answer to this question is relevant 
in the search for a more direct synthesis route. 

A second question concerns the reaction enthalpy of 
||TJ . The ideal hydrogen storage material should produce 
a hydrogen pressure of 0.1 MPa at room temperature. 
The entropy contribution of hydrogen gas at this tem- 
perature favors the right-hand side of Q. At T = 
the hydrogen desorption reaction should therefore have 
an enthalpy of ~ 0.4 eV per desorbed H2 molecule^ Fur- 
thermore, the kinetics of hydrogen adsorption/desorption 
should be sufficiently fast. Finding ways of improving 
the kinetics can start from understanding the bonding in 
Mg(AlH 4 ) 2 , which is determined by the electronic struc- 
ture of the material. 

In this paper we report the results of an ab initio study 
on the properties of magnesium alanate. The structure is 
optimized and the electronic structure is calculated. We 
characterize the bonding in Mg(AlH 4 ) 2 and calculate the 
enthalpy of formation from its elements, as well as the 
enthalpy of the dehydrogenation reaction JIJ . 

Computational methods. Total energies are calcu- 
lated within density functional theory (DFT), using 
the PW91 generalized gradient approximation (GGA) 
functional We use the projector augmented wave 
(PAW) methodi2ii^ and a plane wave basis set, as im- 
plemented in the Vienna Ab initio Simulation Package 
(VASP)>i2*i£ii& The atomic positions and the cell param- 
eters, including the cell volume, are optimized by mini- 
mizing the forces and stresses. A 7 x 7 x 7 Monkhorst- 
Pack k-point mesh is used for sampling the Brillouin 
zoneiii A kinetic energy cutoff of 312 eV is used for the 



lated using a higher kinetic energy cutoff of 700 eV to 
ensure convergence. 

If we calculate reaction enthalpies from total energy 
differences only, we neglect the contributions from atomic 
vibrations. Such contributions are negligible for heavy 
elements, whereas they may be significant for hydrogen. 
For each compound involved in the reaction we calcu- 
late its zero point vibrational energy (ZPVE), ^fi^Z,- ojj, 
resulting from the vibrational modes j in the optimized 
structure. Vibrational frequencies Wj are generated from 
a dynamical matrix, whose matrix elements (i.e., the 
force constants) are calculated using a finite difference 
method^ For the hydrogen molecules we calculate a 
ZPVE of 0.29 eV, in good agreement with the value of 
0.27 eV obtained from the experimental frequency^ We 
also consider the zero point rotational energy (ZPRE) of 
the hydrogen molecules. Assuming that ortho- and para- 
hydrogen are produced in a proportion of three to one, 
the average ZPRE of a hydrogen molecule is 0.011 eV, 
using the energy levels given in Ref. Il9i In summary, the 
reaction enthalpies AH are calculated from 

f 

(2) 

where E^°\ denotes the total electronic energy of the re- 
action products / or reactants i, E^J^ E are the corre- 
sponding ZPVEs, and RE is the ZPRE of all hydro- 
gen molecules involved in the reaction. 

DFT calculations using the common density function- 
als give adequate values for ground state properties such 
as total energies and vibrational frequencies. Excited 
state properties are not given accurately, e.g., the elec- 
tronic band gap is typically underestimated by ~ 50%. 
This in fact stems from an unjustified interpretation of 
the Kohn-Sham eigenvalues of DFT as excitation ener- 
gies. To calculate single-particle excitation energies, one 
should solve a quasi-particle equation using the non-local, 
energy dependent self-energy. The GW technique ap- 
proximates the self-energy by a dynamically screened ex- 
change interaction. Constructing this interaction from 
the orbitals and eigenvalues obtained in a DFT calcu- 
lation with the local density approximation (LDA) is 
called the GoWq approximation. It leads to accurate 
band structures and band gaps for a wide range of semi- 
conductors and insulators i^S GW calculations have also 
been successfully applied to metal hydrides 

We start from an LDA calculation using norm con- 
serving pseudo potentials and a plane wave kinetic en- 
ergy cutoff of 748 eV~ The screened interaction GqWq 
is then calculated using the real space, imaginary time 
formalismiS2iS4 For these calculations we use 150 LDA 
states, a 13 x 13 x 19 real space grid sampling of the unit 
cell, and an interaction cell consisting of 5 x 5 x 4 unit cells. 
The quasi-particle equation is solved while neglecting the 
off-diagonal elements of the self-energy between the LDA 
states. We estimate that the quasi-particle band gap of 
Mg(AlH 4 ) 2 is numerically converged to within ±0.06 eV. 



2 




FIG. 1: (Color online) Crystal structure of Mg(AlH 4 ) 2 ; space 
group P3ml 



TABLE I: Optimized crystal structure of Mg(AlH 4 )o, com- 
pared to the experimental structure at 8 K from Ref. 



Cell parameters 




a 


b 


c 




calc. 


5.23 


5.23 


6.04 




exp. 


5.21 


5.21 


5.84 


Wyckoff positions 




x 


y 


z 


Mg (la) 


calc. 













exp. 











Al (2d) 


calc. 


1/3 


2/3 


0.706 




exp. 


1/3 


2/3 


0.699 


HI (2d) 


calc. 


1/3 


2/3 


0.442 




exp. 


1/3 


2/3 


0.424 


H2 (6i) 


calc. 


0.168 


-0.168 


0.812 




exp. 


0.167 


-0.167 


0.811 
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structure with the Mg atoms on the Cd positions and 
AIH4 tetrahedra on the I positions; the space group is 
P3ml (164) £ The structure basically consists of AIH4 
tetrahedra which form close packed double layers per- 
pendicular to the c-axis, alternated with a layer of Mg 
atoms, as shown in Fig. ^ 

Starting from the experimental structure proposed in 
Ref. |7J, we optimize the atomic positions and cell param- 
eters. As it turns out, for unit cell volumes in the range 
125-150 A 3 the total energy only weakly depends upon 
the volume. We map out the total energy as a function of 
the cell volume. At each volume we optimize the atomic 
positions and the cell shape, and allow for breaking the 
symmetry. Interpolating this energy versus volume curve 
gives a minimum energy at a cell volume of 143.26 A 3 . 
Optimizing the structure at this volume gives the final 
results shown in Table [fl The calculated structure has 
P3ml symmetry and is in good agreement with the re- 
cently obtained experimental structure extracted from 
X-ray and neutron powder diffraction data£ 

The Al-Al nearest neighbor distance within a double 
layer is 3.91 A, whereas the shortest Al-Al distance be- 
tween two double layers is 4.66 A. Mg atoms occupy 
octahedral interstitial sites between two double layers 
with Mg-Al distances of 3.50 A. The A1H 4 tetrahedra 
are slightly distorted, but they retain a threefold rota- 
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tion axis parallel to the c-axis. The Al-Hl and A1-H2 
bond lengths are 1.60 and 1.62 A, and the H1-A1-H2 and 
H2-A1-H2 bond angles are 113.0° and 105.8°. The geom- 
etry of the AIH4 tetrahedra is similar to that found in 
isolated (ATD,) - ions, where the Al and H atoms are co- 
valently bonded. 2 - This geometry is quite different from 
those found in neutral A1H X clusters. 2 ^ The minimum 
Al-H and H-H distances between atoms of different AIH4 
tetrahedra is 3.14 A and 2.63 A respectively, indicating 
that the A1H 4 tetrahedra are clearly separated. The Mg 
atoms are octahedrally coordinated by H2 atoms with 
a Mg-H distance of 1.89 A and H-Mg-H angles of 86.9° 
and 93.1°. This coordination is not unlike that found in 
MgH 2 . 

Electronic structure. Fig. |2Ia) shows the electronic 
density of states (DOS) of magnesium alanate obtained 
from the DFT/GGA calculation. It can be compared 
to the calculated DOS of the lattice of (AIH4)™ tetra- 
hedra, shown in Fig. |2{b) . Here the Mg ions have been 
removed and replaced by a homogeneous positive back- 
ground charge. The similarities between Figs. I^a) and 
(b) demonstrate that the (AlELi) - ions strongly con- 
tribute to both the valence and the conduction bands 
of Mg(AlH 4 ) 2 . Such a dominance of the anion is also 
observed in the simple ionic compound NaCl. 27,28 

Projecting the valence states on atomic orbitals shows 
that Al and H contribute a comparable amount, which 
is a strong indication for covalent bonding within the 
(AlHi) - tetrahedra. The splitting into two valence 
bands, as is most clearly visible in Fig. |2fb), is a rem- 
nant of the splitting between states of s-like (Ai) and 
p-like (T2) symmetry in a single tetrahedron. In an iso- 
lated (A1H 4 )~ ion, the sp-gap is ~ 4 eV. This gap is 
closed to a certain extent by the interaction between the 
(AIH4) - ions, which results in a band dispersion of 2- 
3 eV. It shows that, although the interaction between the 
(AIH4) - tetrahedra is not negligible, it is weaker than the 
interaction within a single tetrahedron. 

If we compare the valence bands of Figs. |2f a) and (b) 
in more detail, we observe a small peak in the DOS of 
Mg(AlH 4 ) 2 , which occurs within the sp-gap mentioned 
above. This peak results from the hybridization between 
H and Mg s states. Hybridization with Mg p and d states 
also gives less clearly visible contributions at higher en- 
ergy. In any case, the H-Mg hybridization is much weaker 
than the H-Al hybridization. 

The DFT/GGA band structure of Mg(AlH 4 ) 2 is shown 
in Fig. [3] It has an indirect band gap of 4.4 eV. The 
bottom of the conduction band is at A and the top of 
the valence band located at maximum of AH at 0.7AH, 
[A = (0, 0, |) , H = (i, i, ±)] . The band gap as ob- 
tained from the GqWq calculation is 6.5 eV. This classifies 
Mg(AlH 4 ) 2 as a large band gap insulator, which is typi- 
cal for ionic compounds. The dispersions of the highest 
valence and the lowest conduction bands in a direction 
along the c-axis are rather small. The direct band gap 
at A is 6.9 eV; the direct band gap at Y is 6.9 eV. More 
generally the dispersion of the GqWq bands is very sim- 
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FIG. 2: (a) Electronic density of states (DOS) of Mg(AlH 4 ) 2 . 
The zero of energy is at the top of the valence band, (b) DOS 
of (A1HJ)2 with a positive homogenous background charge. 




FIG. 3: DFT/GGA band structure of Mg(AlH 4 ) 2 . The zero 
of energy is at the top of the valence band. 

ilar to that of the GGA bands. The total band width, 
6.0 eV, of the GqWo valence bands is ~ 0.1 eV larger 
than the the GGA valence band width. 

The layered structure of Mg(AlH 4 ) 2 , see Fig. ^ does 
not imply that the interactions in the compound are 
strongly anisotropic. The dispersions of the bands in var- 
ious directions are similar, compare e.g. the TA and the 
rM directions [M = (~,0, 0)] . This indicates that the in- 
teractions between the ions within a layer (the ab-planc) 
are comparable to those perpendicular to the layers. 

Formation enthalpy. The formation enthalpy of 
Mg(AlH 4 ) 2 with respect to its elements corresponds to 
the enthalpy of the reaction 

Mg(AlH 4 ) 2 ^Mg + 2Al + 4H 2 ( 5 ). (3) 

Here Mg and Al are in the crystalline phase, whereas 
H2 is in the gaseous phase. For aluminum we use the fee 
structure with a lattice parameter of 4.05 A and for mag- 
nesium we used the hep structure with lattice parameters 
a = 3.21 A, c = 5.21 A. The total energies and ZPE cor- 
rections are given in Table ITT1 From these numbers, the 
reaction enthalpies are then calculated using Eq. J2J). 
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TABLE II: Total energies (with respect to non spin polarized 
model atoms), zero point vibrational energies (ZPVE) and 
zero point rotational energy (ZPRE) in eV/formula unit. 





gTOT 


gZPVE 


gZPRE 


Mg 

Al 

H 2 

MgH 2 
Mg(AlH 4 ) 2 


-1.524 
-3.698 
-6.792 
-8.983 
-36.764 


0.001 
0.004 
0.294 
0.402 
1.520 


0.011 



The reaction enthalpy of © is 0.17 eV per H 2 molecule 
on the basis of total energies. If we include the ZPE, 
the reaction enthalpy decreases to 0.10 eV/H2. Since 
the formation enthalpy is positive, it should, in princi- 
ple, be possible to synthesize Mg(AlH 4 ) 2 from the ele- 
ments. Note that at this energy scale, the contributions 
due to the zero point motions of the hydrogen atoms are 
not negligible. In general, they tend to make a negative 
contribution to the formation enthalpy of the metal hy- 
dride, since the motion of a hydrogen atom in the crystal 
is more confined than in the gas phase. 

To calculate the reaction enthalpy of QJ one also needs 
the optimized structure and ZPVE of MgH 2 . MgH 2 
has space group PA^/mnm (136) and its calculated lat- 
tice parameters are a — 4.51 A and c = 3.01 A. The 
magnesium and hydrogen atoms are at the 2a and 4/ 
(x = 0.304) Wyckoff positions, respectively. As a check 
on the accuracy of these calculations we can also extract 



the formation enthalpy of MgH 2 with respect to its ele- 
ments 

MgH 2 ->Mg + H 2 ( ff ). (4) 

The reaction enthalpy of Q is 0.67 eV/H2 without 
ZPE corrections and 0.57 eV/H.2 with ZPE corrections. 
This value is in reasonable agreement with the value of 

0. 76 eV/H2, which is extracted by extrapolating the ex- 
perimental results to zero temperature^^ 

The calculated reaction enthalpy of is 0.003 eV per 
H2 molecule in the gas phase, without ZPE correction. 
This is negligibly small, but consistent with earlier exper- 
imental dataiiS Moreover, including the ZPE correction 
makes the reaction enthalpy actually slightly negative, 

1. e. — 0.06eV/H2. In any case this number is significantly 
less than the ~ 0.4 eV/H2 which, based on thermody- 
namics, is required to make Mg(AlH 4 ) 2 a good material 
for hydrogen storage. Further investigations are needed 
to see whether e.g. alloying would increase the stability 
of magnesium alanate. 
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